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We present an exact diagonalization study of the half-filled Hubbard model on bipartite quasi- 
one-dimensional lattices. In particular, we emphasize the dependence of the ferrimagnetic ground 
state properties, and its associated magnetic excitations, on the Coulomb repulsion U. 
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I. INTRODUCTION 

In the last years the physics of quasi-one-dimensional compounds has been the object of intense theoretical and 
experimental investigations. In this work we study the ground state (GS) and magnetic excitations of two bipartite 
chains motivated by low-dimensional inorganic [lj and organic ferrimagnets: the AB2 chain in Fig. ^a) and the 
ABC chain in Fig. ^b), respectively. We shall use the simplest approach for interacting electron systems on a lattice 
with N sites, namely the one-band Hubbard Model (HM): 

H=-t ^ c h C i<r + U ^2 n it n il> C 1 ) 

where cj CT (ci a ) is the creation (annihilation) operator for electrons with spin a at site i; t is the hopping parameter 
and U is the on site Coulomb repulsion. For the extremely localized regime (U t) and n = N e /N = 1, where iV e is 
the number of particles, the HM can be mapped onto the Antiferromagnetic Heisenberg Model (AHM) 0: 

w = ^E^-4 (2) 

where J = At 2 /U . We should mention that Macedo et al. Q, using a variety of numerical methods, have derived 
the ferrimagnetic nature of AB2 chains modeled by the HM in the strong and weak coupling regimes. Here we focus 
on the half-filled band case (n = 1) and extend these previous studies to much larger systems using Lanczos exact 
diagonalization technique. 

It is also worth mentioning that the two lattices present three electronic bands in the tight-binding limit (U = 0): one 
flat band at e = for the two chains; and two dispersive ones, e(q) = ±2t^/2cos(q/2) and e(q) = ±t^l + [2 cos(q/2)] 2 , 
with q = 2ttI/(N/3) and I = 0...(A/3) — 1 for the AB 2 and ABC chains, respectively. The flat band is closely 
associated with ferrimagnetic properties of both chains at half-filling 0, [f| . A particular feature of the AB2 chain is 
a local invariance under the exchange of the B sites in any cell / • The eigenvalues of the exchange operator being 
pi = ±1. This symmetry leads to a conserved local spatial parity that in the Heisenberg limit asserts that the B 
sites form either a singlet {pi = 1) or a triplet (pi = — 1) bond state. Therefore for the symmetry sector in which all 
pi = — 1 the spectrum is identical to that of the alternating spin- i /spin- 1 AHM chain 0. 

II. FERRIMAGNETIC GROUND STATE 

A theorem due to Licb and Mattis asserts that the GS of the AHM on a bipartite lattice has total spin 
S g = \Ni — A2I/2, where Aq and N2 are the number of lattice sites at the sub-lattices 1 and 2, respectively. So, if 
\N\ — A2I ~ N the system display unsaturated ferromagnetic (F) GS. The coexistence of both F and antiferromagnetic 
(AF) long-range order implies in ferrimagnetism, as rigorously proved by Tian 0. Another crucial step was provided 
by Lieb 0], w ho p roved that the GS of the HM for n — 1 has S g — |Aq — A2I/2 for any U > 0. In this case, Tian and 
collaborators [T(j further established the ferrimagnetic (FERRI) long-range order of the GS if |iVi — N%\ ~ N. The 
unit cell of the two chains of interest has 2 sites in one sublattice and 1 in the other (see Fig. so that S g = N c /2, 
where N c = N/3 is the number of unit cells. However, although ferrimagnetism is expected, the specific magnetic 
structure of a unit cell may strongly depend on the Coulomb coupling, as well as on quantum fluctuations, and is not 
known a priori from the more general results of the theorems stated above. 
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We can probe the magnetic order through the Magnetic Structure Factor (MSF): 

S(q) = ^e i ^-^{S i -S j ), (3) 
id 

which is related to the zero-field static magnetic susceptibility by x{l = 0) = S(0)/(fcsT), where ks is the Boltzmann 
constant and T is the temperature. The condition for a F (AF) ordered state is that S(0) [S(%)] ~ N, so that in a 
long-range FERRI GS the two conditions must be fulfilled. This is indeed the case for the two chains, as shown in 
Figs. n an d|2 both in the strong and weak coupling limits. 




FIG. 1: Magnetic Structure Factor S(q) for the AB2 (a) and ABC (b) chains in the Heisenberg limit (U ^> t). The size ranges 
from N c = 3 to N c = 10. The insets display the size dependence of the Ferromagnetic [5*(0)] and Antiferromagnetic [S'M] 
peaks. Dashed lines are guides for the eye. 
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FIG. 2: Magnetic Structure Factor S(q) for the AB2 (a) and ABC (b) chains using the HM for U = 2t. The size ranges from 
N c — 2 to N c = 6 for the AB2 chain and to N c — 5 for the ABC. The inset presents the size dependence of the Ferromagnetic 
[S(0)] and Antiferromagnetic [S(7r)] peaks. Dashed lines are guides for the eye. 



Due to the critical nature of both chains at low temperatures, the correlation length £ and x(<? = 0) satisfy power 
law behaviour: £ ~ T~ v and \ ~ T" 1 as T — > 0. Since ^ ~ N at T = 0, using scaling arguments and the results 
of Fig. ^ we have T f ~ r -t/ /T, i. e., *y — 1/ = 1, in agreement with the values 7 = 2 and z/ = 1 derived using 
renormalization group techniques pd| . 

In Fig. |31we present the local magnetization (5f) and the average local value of double occupancy ((niiriii)) as 
function of U. Since the system has a spontaneously broken SU (2) symmetry, we can choose the spin sector S z = S g 
to study these quantities, where S z is the total spin z component. In the Heisenberg limit (U 3> t, (Tii-|-7ij^)=0), the 
bulk magnetization estimated from chains up to N c = 10 are: (S\) = —0.2925 and (Sj^ ) = 0.3962, for the AB2 
chain, in agreement with density matrix renormalization group (DMRG) calculations [12| and with values for (S\) 
and 2(S z Bi 2 ) in the spin-i/spin-1 chain [l^. For the ABC chain, (S%) = -0.2567, (S Z B ) = 0.3741 and (Sg.) = 0.3826. 



3 



Notice that a FERRI structure sustain quantum fluctuations, which strongly reduce the local magnetization from the 
classical value, although in any case the unit cell magnetization remains 1/2. In the non-interacting limit (U = 0) 



we see in Figs. Efa) and|3fb) that {S Z A ) = and (S 



Bi(B) 



0.5 for the AB 2 {ABC) chain. These results 



filling the bands up to n = 1 and choosing the spin sector S z — S g , there exist 



being related to the band structure: nmug uk uanuo up uu ID - i an^ uiuuomg one ocuvi u — u g , 
N c unpaired electrons in the flat band which occupy B\(B) and ^(C 1 ) sites only. It is important to note that in this 
limit the global GS does not have an unique value of S, but rather S < N c /2. As the Coulomb coupling is turned on, 
the GS, with a unique S g = N c /2, evolves adiabatically from the weak to the strong coupling regime (Fig. [3Jl. 
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FIG. 3: Site magnetization (Sf) and average site double occupancy (nifTij^) as function of U for N c 
chains in the spin sector S z = S g . Dashed lines are guides for the eye. 
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III. MAGNETIC EXCITATIONS 



Systems with a FERRI GS naturally have F (lowering the GS spin) and AF (rising the GS spin) magnons as 
their elementary magnetic excitations. The two chains have three spin wave branches: an AF mode, defined as 
AE s+ (q) = E(S Z = S g + l,q)-E ; and two F ones, derived from AE S -{q) = E(S Z = S g -l,q)-E , where E is the 
GS energy and E(S z ,q) are lowest excitation energies in the sector {S z ,q}, with q = 2nl/N c and I — 0, 1, N c /2. 
The acoustical and optical F modes are identified from symmetry considerations. These modes arc depicted in Fig. 
IIIII for the two chains using the AHM: the AF mode has a gap, denoted by A s+ ; the Fl mode is gapless, i. e., the 
Goldstone mode, consistent with the symmetry broken phase of the chains; the gap of the F2 mode is denoted by 
As- 

For the AB 2 chain the gapped F2 branch is flat, Fig. Wlf al . and the other modes are dispersive. Since these 
dispersive modes preserve the local triplet bond, they are identical to those found in the spin-^/spin-l ^| chain. 
Surprisingly, Linear Spin Wave Theory (LSWT) [1J| predicts that As- = I, very close to our estimated value: 
As- = 1.0004J. Moreover, a good agreement is found for the gapless Fl branch in the long wave-length limit. 
However, both LSWT and mean field (MF) theory jjjjjj predicts As+ = 1, deviating from our estimated exact 
diagonalization value: Ag + = 1.7591J, which is in excelent agreement with numerical and analytical calculations for 
the spin-i/spin-1 chain ^S]. On the other hand, the Interacting Spin Wave Theory 0] derives a better result for 
As + , but it implies in a higher shift for As- (flat mode) not observed in our data of Fig. Illlf a) . This flat mode 
indicates the localized nature of the excitation, which is associated with the formation of a singlet state (even parity 
symmetry) between the B sites in one cell, while the other cells have odd parity symmetry (B sites in a triplet state) . 
Thus, since the symmetries of the F modes (dispersive and flat) are distinct, a possible level crossing is not avoided, 
as in fact observed in the data of Fig. Illlf ah for q > 0.757T the localized excitation is below the dispersive one. 

For the ABC chain we can see in Fig. IHIf b) that the three branches are dispersive. Like in the AB 2 chain, the 
AF mode is gapped with A s+ = 1.4189. The LSWT [l7| conducts to a qualitative but not quantitative correct AF 
branch. However, a good agreement is found for the gapless Fl mode in the long-wavelength limit. Moreover, the gap 
for the F2 mode is estimated to be As- — 0.6778 J. We can clearly observe in Fig. Illlf b) the level repulsion between 
the F branches, which in this case have the same symmetry. 
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FIG. 4: F and AF spin wave modes of the AB2 (a) and ABC (b) chains for N c = 10 (circles), 8 (triangles down), 6 (triangles 
up). Solid lines are the LSWT results from Ref. [3 (AB 3 ) and Ref. [13 (ABC); dashed lines are guides for the eye. 



IV. MAGNETIZATION PLATEAUS AND SPIN WAVE GAPS 

The AF gaps found above are responsible for plateaus in the magnetization per spin [m(H) — (S z ) /N] as function 
of the applied magnetic field H . In fact, it has been shown |l8j | that if f(s — m) = integer, a plateau may appear in 
the magnetization curve of the AHM. In the last equation, s is the site total spin, to is the magnetization and / is the 
number of sites in one unit period of the GS for a given value of H . The Hamiltonians that we are considering here 
have s = 1/2 and three sites per unit cell, so, unless the system spontaneously breaks the translation invariance (e.g., 
by growing the unit cell), we expect a plateau at m = 1/6. This is indeed observed in Fig. [31 and its width matches 
exactly As+, which measures the region of stability of the FERRI phase. For higher fields the magnetization increases 
in the expected way [l^, as shown in the inset (2) of Figs. El (a) and [51(b), before the magnetization saturates at 
to = 1/2 for H = 3 J and H = 2. 3660 J, respectively. This field dependent behaviour contrasts with the linear one 
predicted by MF theory [l5j. 
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FIG. 5: Magnetization as function of the applied magnetic field H for the AB2 (a) and ABC (b) chains using the AHM. The 
inset (1) presents a finite size scaling study for m = 1/3; the inset (2) is a curve traced from the midpoints of the steps using 
the finite size results, except at plateau regions. 



Now we turn to study the evolution of the AF and the F2 spin wave gaps as function of U using the HM: Ag + (?7) 
(plateau width) and As-{U), respectively. In order to reduce finite size effects in calculating As±(U) we have used 
the boundary condition that minimizes the gap between the GS and excited levels. Our estimated curves are depicted 
in Fig. Housing the scaling ansatz 0] for extrapolation: A = Aoo + ciN~ 2 + C2N~ 4 , with N — 9, 12, 15 and 18 for the 
AB 2 chain; N = 6, 9, 12 and 15 for the ABC chain. Note that A S +{U = 0) = and 2t for the AB 2 and ABC chains, 
respectively. This can be understood by noticing that As+(U = 0) is given by the gap between the two dispersive 
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tight-binding electronic bands. Moreover, As+(U — > 0) ~ U x , x ~ 2, for the chain, whereas for the ABC chain 
Ag + (?7 — > 0) increases linearly before reaching the tight-binding limit 2t. For both chains, we observe a crossover 
to the Heisenberg limit: ~ J, in agreement with the results of Fig. IIIII Further, we notice that As~(U — » 0) 

nullifies linearly with U, for both chains, following the behaviour of the gap between the flat electronic bands 0. 
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FIG. 6: Extrapolated AF (As+) and F2 (Ag_) spin wave gaps in units of t for the AB2 (a) and ABC (b) chains as function 
of U using the HM. Solid (dashed) line is the AF (F2) gap for the AHM with J = At 2 /U; dotted lines are guides for the eye. 



V. CONCLUSIONS 



In this work we reported a detailed exact diagonalization study of the half-filled Hubbard model on bipartite 
quasi-one-dimensional lattices. These chains have experimental motivation in inorganic and organic ferrimagnetic 
polymers. We have focused on the description of their ferrimagnetic ground state properties and magnetic excitations. 
In particular, we have revealed the non-trivial physical behaviour exhibited by these chains as function of the Coulomb 
coupling strength. 
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